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Differential magneto-optics are used to study the effect of dc in-plane magnetic field on hysteretic 
behavior due to geometrical barriers in Bi2Sr2CaCu208+5 crystals. In the absence of an in-plane field 
a vortex dome is visualized in the sample center surrounded by barrier- dominated flux- free regions. 
With an in-plane field, stacks of Josephson vortices form vortex chains which are surprisingly found 
to protrude out of the dome and into the vortex- free regions. The chains are imaged to extend up to 
the sample edges, thus providing easy channels for vortex entry and for drain of the dome through 
geometrical barrier, suppressing the magnetic hysteresis. Reduction of the vortex energy due to 
crossing with Josephson vortices is evaluated to be about two orders of magnitude too small to 
account for the formation of the protruding chains. We present a model and numerical calculations 
that qualitatively describe the observed phenomena by taking into account the demagnetization 
effects in which flux expulsion from the pristine regions results in vortex focusing and in the chain 
protrusion. Comparative measurements on a sample with narrow etched grooves provide further 
support to the proposed model. 

PACS numbers: 74.25.Ha, 74.25.0p, 74.25.Uv, 74.25.Wx 



I. INTRODUCTION 

One of the most prominent characteristics of type II su- 
perconductors is the large magnetic hysteresis that arises 
from the existence of potential barriers that impede vor- 
tex motion.^ There are several sources of potential barri- 
ers: pinning due to quenched material disorder in the 
bulk of the sample, surface pinning produced by the 
roughness of the sample surface,^ Bean-Livingston sur- 
face barrier arising from the force between a vortex par- 
allel to a surface and its image, ^ and geometrical barrier 
(GB) that arises from the force due to Meissner currents 
in thin slab-shaped superconducting samples in perpen- 
dicular field^^ 

Strong vortex pinning and the associated magnetic 
hysteresis are vital elements in superconductor applica- 
tions that are based on high critical currents. The mag- 
netic hysteresis, however, masks the equilibrium magne- 
tization properties, the study of which is essential for 
comprehension of the thermodynamic structure of the 
vortex matter. In order to overcome this limitation, vor- 
tex shaking has been employed in recent years, which was 
shown to suppress the magnetic hysteresis very effectively 
at not too low temperatures. ^^^^ In this method an ac 
in-plane field is applied in addition to the main dc 
Hz which causes local shaking and equilibration of the 
vortices, dc in-plane fields were also shown to reduce the 
magnetic hysteresis^^^-iS. 

Vortex shaking was found to suppress the hysteresis 
that arises either from bulk pinning that is usually the 
dominant source of hysteresis at lower temperatures, sur- 



face barriers that often dominate vortex dynamics at el- 
evated temperatures and fields, ^^^^ and GB's that com- 
monly govern the hysteresis at high temperatures and low 
fields J^i^ii^i2ii2^ Detailed theoretical studies explained the 
suppression of bulk pinning hysteresis by vortex shak- 
ing in terms of the Bean model in a turning magnetic 
field. ^^^^ This description may similarly apply to the 
case of surface pinning. Also the suppression of the 
Bean-Livingston surface barrier hysteresis could possibly 
be understood in terms of easier vortex loop nucleation 
at the sample corners for tilted vortices. In the above 
mechanisms the potential barriers are microscopic and 
hence local ac displacement or tilting of a vortex seg- 
ment in the presence of an average dc driving force can 
enhance vortex activation and relaxation. In contrast, 
the GB in a platelet sample gives rise to a vortex dome 
in the center of the sample surrounded by a vortex-free 
region. This vortex-free region constitutes a large bar- 
rier that has a macroscopic energy scale of the order of 
the vortex line tension times the sample thickness and a 
macroscopic length scale of the order of sample's width. 
As a result, this barrier cannot be overcome by any local 
vortex perturbation or tilting. Hence the experimental 
observation that the vortex shaking by the in-plane field 
suppresses the hysteresis due to GB was highly surprising 
and has remained unresolved so far. 

In this paper we describe a possible resolution of this 
question in the case of layered superconductors by in- 
vestigating the GB behavior in the presence of a dc in- 
plane field. In highly anisotropic layered materials like 
Bi2Sr2CaCu208+5 (BSCCO) the in-plane field H^, re- 
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suits in the formation of stacks of Josephson vortices 
(JV's).^^^^ The addition of causes the formation 
of crossing lattices in which stacks of pancake vortices 
(PV's) form chains along the Josephson vortices^^ — 
Surprisingly, we find that these chains exist also out- 
side the vortex dome, protruding like vortex whiskers 
into regions that are vortex free in the absence of H^. 
As a result, easy channels are apparently formed that 
allow vortex penetration through the macroscopic GB 
suppressing the hysteresis. This observation appears to 
be inconsistent with theoretical expectations since the re- 
duction in the energy of the pancake stack by the Joseph- 
son vortices should be much smaller than the height of 
the geometrical barrier. ^'^^ Therefore, the regions out- 
side the vortex dome should have remained vortex free 
even in the presence of and no significant suppression 
of the barrier should occur. We propose an explanation 
to this phenomenon in terms of interaction between JV's 
and PV's that is enhanced by the demagnetization ef- 
fects in platelet geometry. We demonstrate this scenario 
by studying vortex behavior in a sample with a shallow 
and narrow groove etched into its surface using a focused 
ion beam (FIB). 



II. GEOMETRICAL BARRIER 

The GB is a macroscopic potential barrier that vortices 
have to overcome in order to enter or leave the sample 
and is the main source of magnetic hysteresis in thin slab- 
shaped samples with weak pinning. The barrier is formed 
by the interplay between the vortex line tension and the 
Lorentz force that is induced by the circulating Meissner 
currents. In a sample with elliptical cross section, for 
example, the energy of a test vortex has two contribu- 
tions: the positive vortex line energy increases gradually 
from zero at the edge of the sample to a value of eod in 
the center (eo is the vortex line energy per unit length 
and d is the sample thickness). The Meissner currents 
that fiow over the entire surface exert inward forces on 
the vortex that gradually lower its energy as it moves to- 
ward the center. In a sample with elliptical cross section 
the two contributions cancel out rendering a constant 
vortex energy and a corresponding position independent 
and reversible induction B{x^ H). In a platelet sample of 
rectangular cross section, in contrast, by cutting through 
the sharp rims, the vortex attains its full line energy eod 
within a distance of the order of d/2 from the edge,^ while 
the Meissner currents remain distributed as in elliptical 
sample. As a result, a geometry related barrier of height 
of the order of eod is formed that extends over a width 
of the order of the half width of the sample, wi^ —i^ 
Consequently, when the field Hz i s inc reased above the 
penetration field Hp {2Hci/'n)^ d/w vortices entering 
through the edges are swept by the Meissner currents to- 
ward the center where they accumulate giving rise to a 
dome shaped induction profile of width 26, see Fig. [IJc). 
The vortex-filled dome is current free while the surround- 



ing vortex- free region carries the Meissner currents. The 
induction profile, with a central dome surrounded by an 
area of zero B is very distinctive and is the hallmark of 
the GB model. The vortices in the dome cannot leave 
the sample since they are trapped by the barriers due 
to the Meissner currents in the surrounding fiux-free re- 
gions. Thus, when the field is decreased the dome ex- 
pands, while keeping the total fiux in it constant, giving 
rise to hysteretic magnetization. Only when the dome 
reaches the edges the barriers vanish allowing vortices to 
leave the sample. 



III. EXPERIMENTAL DETAILS 

Several over-doped and optimally-doped BSCCO sam- 
ples were studied. We present here data of two over- 
doped samples, with a critical temperature of Tc = 87K, 
cleaved from one thick crystal, sample A and sample 
B, with platelet shapes of length 2L = 1500/im, width 
2w = 420/im and thicknesses d of 8 and 20/im re- 
spectively. Similar results were obtained from opti- 
mally doped crystals. The differential magneto-optical 
(DM0) technique^^ was used to study the induction pro- 
files and the magnetization behavior of the samples. By 
minimizing the distance between the magneto-optical in- 
dicator and the sample, and by utilizing differential imag- 
ing we achieved high visibility of the vortex chains and 
studied their effect on hysteretic magnetization. The to- 
tal exposure time of a differential image was typically 2 
to 3 minutes. We found that following a field ramp the 
JV's often show slow relaxation for about two minutes, 
during which their position within the sample changed. 
Thus quality images of the chains could be achieved by 
waiting 2 to 3 minutes between the end of a field ramp 
and the beginning of image exposure. In order to sim- 
ulate the effect of JV's, we also studied the magnetic 
behavior of crystal A after two fine grooves, 80 and 160 
nm deep and 350 nm wide, were etched on the surface 
using an FEI Helios Focused Ion Beam system. 



IV. RESULTS 

In order to visualize the vortex dome structure with en- 
hanced sensitivity, we carried out DM0 with either field 
modulation or current modulation. In the former case, 
Bix^y) images taken at Hz = Ho — dH are subtracted 
from images taken at Hz = Ho-\- dH and few tens of such 
differential images are averaged resulting in the differen- 
tial image SB{x^y).^^ When the field is modulated by 
2dH the vortex dome expands and shrinks periodically, 
keeping the total fiux in the dome constant; see Fig. [U^c). 
As a result, when the field is raised to Hz = Ho ^ dH 
the induction in the center of the dome is increased giv- 
ing rise to a positive 5B signal, while the induction near 
the dome edges is decreased resulting in a negative 5B 
there, see Fig. [TJd). This effect is clearly visible in 



3 



2D 


i ^ 

— H^=Hp+dH 
















n- 






■1 -O.S 

1 


0.5 1 




-200 -tKJ 100 2W 








Figure 1. (color online) (a) Field modulated DM0 image of 
sample B displaying the GB dome along the central part of the 
sample {H, = 6.7 Oe, dH = 0.4 Oe, H^: = 0, and T = 75 K). 
The dark (white) vertical stripes indicate the dome (sample) 
edges, (b) The differential field profile SB(x) across the width 
of the sample averaged over the strip marked by the dashed 
lines in (a). The two minima with negative 6B indicate the 
dome edges, and the two maxima - the sample edges, (c) The 
theoretical field profiles of a thin and infinitely long sample 
B{x,z = 0), at two values of = Hq ± dH (Ho = 6.7 
Oe, dH = 0.4 Oe, Hd = 33.1 Oe, d/w = 0.1). Inset: The 
experimental schematics showing a BSCCO crystal with gold 
contacts glued onto a substrate, (d) The difference 6B{x) 
of the two B{x) profiles plotted in (c) corresponding to the 
experimentally measured differential signal. The dome width 
26 is marked. 



the DM0 data in Figs. [11(a) and (b), in which the edges 
of the dome appear as dark regions of negative values. 

In DM0 with current modulation, images taken in the 
presence of transport current I = —dl are subtracted 
from images taken at / = and averaged. Since the 
Meissner currents flow in the opposite directions on the 
two sides of the dome, the applied transport current in- 
creases the total current flowing on one side of the dome 
and reduces the current on the other side. As a result, 
the dome is shifted away from the center as shown in Fig. 
[2jc). Therefore, in the current modulated DMO images 
one edge of the dome will appear as a bright strip with 
a positive signal and the other as a dark strip with a 
negative signal as shown theoretically in Fig. [2](d) and 
experimentally in Figs. [2] (a) and (b). Note also that the 
self-induced magnetic field of the transport current has 
an opposite sign at the two edges of the sample. As a 
result in DMO image with current modulation one edge 
of the sample is dark while the other is bright. Fig. [Ifa), 
in contrast to DMO with field modulation in which both 
edges appear bright. Fig. [IJa). 

We now analyze the behavior of the dome width, 26, 
along the magnetization loop. Figure [31(a) shows the cal- 
culated behavior of b using the GB formalism of Refs. 
0, HSi, and Initially the sample is in the Meiss- 



Figure 2. (color online) (a) DMO image with current modu- 
lation of sample B displaying the GB dome along the sample 
center at H^ 6.8 Oe, if^ = 0, T = 75 K, and dl 10 
mA. The current modulation causes the right (left) edge of 
the dome to appear bright (dark) and the right (left) edge of 
the sample to appear dark (bright), (b) The differential field 
profile across the width of the sample averaged over the strip 
marked by the dashed lines in (a), (c) Theoretical field pro- 
files across an infinitely long and thin sample in the presence 
of positive and negative applied transport current showing 
the shift of the dome, (d) The difference of the two profiles in 
(c) corresponding to the experimentally measured differential 
signal. 



ner state (b = 0) which is retained until Hz = Hp 
{2Hci/'n)\/d/w^ above which flux starts entering the 
sample causing b to increase rapidly. As Hz is further in- 
creased, b asymptotically approaches its maximum value 
^max — w — corresponding to the dome width at 
which the potential barrier vanishes. If the field is de- 
creased from an intermediate value, b first expands to 
^max while the total flux in the dome remains constant 
(open circles in Fig. [3](a)). When the field is further re- 
duced, b remains at its maximal value allowing the vor- 
tices to exit through the edges. The corresponding cal- 
culated local magnetization loop is shown in Fig. [D^a). 

It can be readily shown^ that both the ascending 
and descending magnetization branches are energetically 
metastable. Compared to the equilibrium conditions, at 
which the vortices in the dome have the same energy 
as the vacuum £^ = 0, on increasing (decreasing) field 
the vortex energy is negative (positive), the dome is nar- 
rower (wider), there is a deficiency (excess) of vortices in 
the sample and the magnetization is lower (higher) than 
at equilibrium. The dashed lines in Figs. [3{a) and [U^a) 
show the corresponding equilibrium behavior of the dome 
width and of the local magnetization. 

Figure [Sl^b) shows the experimentally measured dome 
width along the magnetization loop in sample B. In the 
absence of an in-plane field {H^ = 0) the behavior is qual- 
itatively very similar to the theoretical GB curve shown 
in Fig. [31(a). On increasing Hz the sample is initially in 
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the Meissner state {b = 0) until > Hp = 4.5 Oe where 
b starts increasing rapidly. On decreasing from its 
maximum value (10 Oe), b continues to increase initially 
until it reaches a maximum value, which is then roughly 
retained on further field decrease. This hysteretic be- 
havior of the GB is markedly altered by the presence of 
even a small in-plane field. For = 5 Oe the pene- 
tration field Hp is slightly reduced to 4.1 Oe and the as- 
cending field branch is similar to the H^ — case. How- 
ever, on the descending branch, b is significantly below its 
maximum value and decreases rapidly, reaching zero at 
finite H^. The shrinkage of the dome width on decreas- 
ing field is in sharp contrast with GB model according 
to which vortices are trapped in the dome as long as the 
surrounding vortex- free region is present {b < w — 
and can leave the sample only when dome reaches the 
sample edges {b = w — d/2). This indicates that in the 
presence of in-plane field the vortices can surprisingly exit 
the sample by permeating through the vortex-free region 
that is impenetrable within the GB model. At a higher 
in-plane field Hx = 30 Oe the dome is formed already at 
Hz = 2 Oe, it shrinks back to zero with decreasing field 
for Hz < 2 Oe, and the dome width b behaves essen- 
tially reversibly following the equilibrium curve (dashed) 
of Fig. [3Ka). This shows that the vortices in the dome 
are effectively in equilibrium with the vacuum and unex- 
pectedly can now readily move in and out of the dome 
through the vortex-free region as Hz is varied. 

The DM0 imaging allows a very sensitive measure- 
ment of the dome width, but cannot provide a measure 
of the dc magnetization. In order to analyze the sam- 
ple magnetization we have carried out a dc MO study in 
sample A. Due to the lower sensitivity the resulting data 
in Fig. m^b) are nosier and less accurate, but qualita- 
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Figure 3. (color online) (a) The calculated half- width of the 
vortex dome b as Hz is increased from zero to Hz ^ Hci and 
decreased to zero again (solid line) within the GB model for 
an infinitely long strip with d/w = 0.1. When Hz decreases 
from an intermediate value, the dome widens following a con- 
stant flux path (open circles) until it reaches the maximum 
value oi b — w — d/2 at which the flux can exit the sample. 
The equilibrium dome width (dashed line) is shown for com- 
parison describing the conditions of zero vortex energy in the 
dome, (b) The dome width measured in sample B at T = 75 
K while ramping Hz from zero to 10 Oe and back in the pres- 
ence of zero and non zero Hx showing that the in-plane fleld 
equilibrates the GB hysteresis. 
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Figure 4. (color online) (a) The calculated local magnetiza- 
tion loop B — H in the sample center within the GB model for 
the same parameters as in Fig. [3l^a) (solid line) as compared 
to equilibrium behavior (dashed). AB is the width of the 
hysteresis, (b) The measured B — H in the center of sample 
A (before FIB etching) at H^ = 0, 40 Oe and T = 79 K. 



tively consistent. Figure [3Ka) shows the calculated local 
magnetization loop B — H in the center of the sample 
(solid lines) along with the equilibrium magnetization 
(dashed line). The large hysteresis is the result of the 
GB that gives rise to wide fiux free regions on ascending 
field and essentially no fiux-free regions on descending 
field. The experimentally measured local magnetization 
loop using MO imaging without modulation is shown in 
Fig. [11(b). In the absence of an in-plane field {H^ = 0) 
the behavior is consistent with the GB model. Applica- 
tion of Hx^ however, significantly reduces the hysteresis 
on the ascending and descending branches showing that 
Hx facilitates both vortex entry and exit. 

In order to quantify the reduction of the magnetic hys- 
teresis we plot in Fig. [5] the measured hysteresis width, 
\AB\ (see Fig. HJa)), as a function of Hz and Hx- The 
figure shows that as \Hx\ is increased the hysteresis is 
suppressed significantly, however, the loop does not close 
completely and some hysteresis is seen around the pen- 
etration field, even at the highest in-plane field we mea- 
sured, 140 Oe. 

The evidence shown so far points to the fact that the 
JV's introduced by the in-plane field significantly reduce 
the effectiveness of the geometrical barriers for some rea- 
son. In order to comprehend the microscopic origin of 
this phenomenon we have optimized our DM0 imaging 
technique to resolve the PV chains along the JV's. The 
stacks of PV's that decorate the stacks of JV's were pre- 
viously observed by Bitter decoration, ^^^^^ scanning Hall 
probes^^'^^ and MO imaging. ^^^^'^^^^ We find that by ad- 
justing the modulation parameters the visibility of the 
chains can be greatly enhanced using either field or cur- 
rent modulation. This is demonstrated in Fig. [6l where 
the vortex chains are clearly visible as bright slightly 
tilted horizontal lines that are aligned along the direc- 
tion of the in-plane field Hx. We verified that these lines 
are indeed vortex chains decorating the stacks of JV's 
by rotating the in-plane field and by varying Hx. We 
find that the lines are always aligned along the direction 

of Hx and their separation follows \/2j^o/\^\Hx\ with 
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= 480 ± 70, where = 20.7 G/im^ is the flux quan- 
tum and 7 is the anisotropy parameter. In Fig. [6] the 
vortex dome is visible as a half-oval envelope that has a 
bright contrast on the right-hand-side and dark on the 
left-hand-side due to the current modulated DM0. We 
expect PV's to be present only inside the dome, residing 
on the JV stacks and between them depending on fleld 
strength, while in the vortex-free regions outside the 
dome the PV's should be absent due to the large Meissner 
currents and correspondingly high vortex potential!^ — 
The striking observation in Fig. [6] is that the vortex 
chains are visible both inside the dome as well as outside 
the dome forming whisker-like protrusions that extend 
all the way to the sample edges. This means that in the 
presence of an in-plane fleld PV's can reside in the regions 
outside the dome that are vortex free in the absence of 
Hx. A similar situation was observed by decoration ex- 
periments on BSCCO crystals in Ref. J^. In Fig. [6] the 
vortex chains inside the dome have a stronger contrast as 
compared to the chains outside the dome. We flnd that 
this is not always the case and the whiskers may have 
stronger contrast than the chains in the dome depending 
on the fleld and modulation parameters. 

We now present the hysteretic behavior of the dome 
with and without Hx- Figure [71(a) shows a DM0 image 
using current modulation at i^^ = 3 Oe in the absence of 
an in-plane fleld. At this fleld the sample is in the Meiss- 
ner state on the ascending branch and the vortex dome 
is absent, since Hz < Hp. On the descending branch, 
however, the sample is in the mixed state and the dome 
is maximal as shown in Fig. [71(b). This highly hysteretic 
behavior of the dome is the result of flux trapping due to 
GB. In contrast, in the presence of H^ = 30 Oe the dome 
has the same width on the ascending and the descending 
flelds, showing that the vortices are in equilibrium state. 
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Figure 5. (color online) The width of the magnetization loop, 
|AB| (see Fig. [J}, measured by sweeping Hz at various values 
of Hx' The discontinuities about Hx ~ are an artifact due 
to the MO indicator. 
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Figure 6. A current modulated DMO image of the decorated 
JV's (slightly tilted horizontal lines) inside and outside the 
dome in sample B. The two side edges of the sample are 
marked by gray arrows and the dome edges by black arrows. 
The center of the sample coincides approximately with the im- 
age bottom while the top edge of the sample is slightly above 
the top of the image. As a result, the visible dome has the 
shape of half an ellipse. The decorated JV's display a large 
contrast inside the dome while their whisker- like protrusions 
outside the dome display a lower contrast. The image was 
taken at T = 77 K, Hz = 3.6 Oe, dl — Q mA and the gray 
scale corresponds to about 0.5 G. The image width is 0.4 mm. 



Figure [8] shows similar DMO images using fleld modula- 
tion at Hz = 5 Oe. In absence of an in-plane fleld Fig. 
[81(a) shows a narrow dome just above the penetration on 
ascending fleld, while Fig. [Sfb) shows a very wide dome 
on the descending Hz branch. When H^ = 30 Oe is ap- 
plied, the domes on ascending and descending branches 
become essentially identical. A movie showing the com- 
plete magnetization loop from which the Fig. [8] is ex- 
tracted is available in Ref. |47|. 

The energy of the PV's that reside on JV's is lower due 
to their attractive interaction. Since the JV's transverse 
the entire sample, one could argue that the energy of the 
PV's could be lowered sufficiently by the JV stacks to 
extend into the GB in the vortex free region. The GB 
energy scale is Eg^ = eo<i, where eo = i^ci^o/47r, with 
the barrier being highest near the sample edges and grad- 
ually decreasing toward the dome.^ The JV-PV interac- 
tion strength per unit length of the PV stack, is given by 
crossmg energy^ 

^""^V 2$o 47rVsln(3.57s/A)' ^' 

where 7 = 480 as described above, s = 15 A. is the layer 
spacing, and A is the London penetration depth such that 
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Figure 7. Current modulated DM0 images at Hz — 3 Oe 
{dl — 10 mA, T = 75 K) of sample B taken along the mag- 
netization loops shown in Fig[3l^b) for Hx = (left column) 
and Hx — 30 Oe (right column) on ascending (top row) and 
descending (bottom row) Hz branches. The left edge of the 
dome is dark and the right is bright (marked by black arrows) , 
while the left edge of the sample is bright and the right is dark 
(marked by gray arrows). The slightly tilted horizontal stripes 
in (c) and (d) are PV chains. The gray scale is about 2.5 G 
in (a) and (b) and 1.5 G in (c) and (d). The images are 0.4 
mm wide. 
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Figure 8. Field modulated DM0 images at if^ = 5 Oe {dH — 
0.4 Oe, T = 75 K) of sample B taken along magnetization 
loops in the presence of Hx = (left column) and Hx = 30 
Oe (right column) on ascending (top row) and descending 
(bottom row) Hz branches. The dome edges are dark (marked 
by black arrows) and sample edges are bright (marked by gray 
arrows). The slightly tilted horizontal stripes in (c) and (d) 
are PV chains. The gray scale is about 1.25 G, 2.5 G, 1.8 G 
and 1.9 G in (a) to (d) respectively. The images are 0.5 mm 
wide. A movie showing the complete magnetization loops is 
available in Ref. 47. 



A = A(0)/v'l - {T/TcY with A(0) = 2000A. In the case 
of our experimental conditions the resulting total vortex 
energy reduction is e^d 10~^£^^5. Therefore, the de- 
gree of the spacial extension of the PV whiskers into the 
vortex-free region should be extremely small, O(10~^) of 
the width of the flux-free region, and hence should be 
hardy observable. Nevertheless, our images clearly show 
that PV chains extend through the entire vortex-free re- 
gion, in sharp contrast to the above estimates. 

It was recently shown that small local inhomogeneities 
in Hci could lead to very large variations in the equilib- 
rium local induction, a phenomenon attributed to demag- 
netization effect s.^^ In order to see whether demagneti- 
zation can explain the phenomenon at hand we studied 
the GB mechanism in sample A after two shallow grooves 
were etched on its surface across the width of the sample 
as described above. We found that the grooves have a 
surprisingly large effect: a 160 nm deep groove, just 2% 
of the sample thickness, reduces the GB hysteresis com- 
parable to the suppression by the in-plane field of about 
10 Oe, and vortex whiskers are formed along the groove. 
The current modulated DMO image in Fig. [9] shows that 
on ascending field the whiskers along the groove (marked 
by arrows) have a contrast that is comparable to the 
dome in the pristine region but their extent is about 60% 



wider than the pristine dome. On the descending Hz 
branch the vortex dome in the pristine region closes grad- 
ually while the whiskers along the groove extend up to 
the sample edges similar to the behavior of the whiskers 
along the JV's. 



V. DISCUSSION 

The reduction of the GB by JV stacks has four main 
manifestations in our data: {%) suppression of the mag- 
netic hysteresis, (n) reduction of the penetration field 
with increasing \Hx\^ (in) the presence of PV chains out- 
side the bulk vortex dome, and (iv) on the descending 
Hz branch of the magnetization loop — the reduction of 
the bulk dome width to zero at some finite Hz = Hq, 
with Hq increasing with increasing \Hx\. A shallow and 
narrow groove in the absence of H^, had the same qual- 
itative effects on the GB. Indeed, both the groove and 
a stack of JVs reduce the energy of the vortices residing 
on them: A groove reduces the average vortex energy by 
eod/d, where S is the groove depth, while a stack of JVs 
reduces the vortex energy by ex- Thus in both cases a 
chain-like structure with enhanced density of vortices is 
expected to be formed along the groove or the JVs. In 
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Figure 9. (color online) A zoomed-in view of the vortex dome 
in the central part of sample A after a groove 160 nm deep 
and 350 nm wide was etched across the width of the sample 
using FIB (marked by arrows). The image was obtained using 
current modulated DM0 on ascending Hz = 4 Oe at T = 82 K 
and dl = 10 mA. The dome left (right) edge is dark (bright) 
and its schematic induction profile is shown by the dot-dashed 
line. Vortex whiskers protruding outside the dome are visible 
along the groove. 



the following discussion, we therefore use the term groove 
in order to describe the related reduction of the geomet- 
rical barrier in both experimental situations. Note that 
in contrast to the fixed groove fabricated by FIB, in the 
case of JVs both the density of the "grooves" and their 
effective "depth" de^/eQ increase with (see Eq. ([T])). 
As a result, in the case of JV "grooves" one may expect 
a stronger suppression of the GB with increasing as 
observed experimentally. 

In the discussion that follows we show that the sup- 
pression of the vortex energy along the groove is insuffi- 
cient to account for the observed phenomena on its own. 
We then demonstrate that the boundary conditions re- 
quire the presence of substantial shielding currents that 
flow along the groove. These currents cause an additional 
significant increase of the magnetic induction inside the 
groove enhancing the suppression of the GB in accor- 
dance with the experimental results. 

We start by evaluating the expected reduction of the 
GB due to a groove across the width of the sample. 
Let us consider a sample of width 2w and thickness d 
{—d/2 < z < d/2) having a groove of depth 5 extend- 
ing from X = —w to X = w. Since the vortex energy 
in the groove is lower, we expect that the vortex dome 



along the groove will be wider than in the bulk forming 
whisker-like protrusions of vortices into the flux free re- 
gion. If the groove width, denoted by 2/, is large enough, 
i.e. / > w then it can be described by the GB formu- 
lation following for example Ref. 43. Figure [10] shows 
for example the resulting vortex potential U {x) and field 
B{x) profiles in the bulk and along the groove of depth 

5 = 0.1 3(i at Hz = 0.065i^ci under equilibrium condi- 
tions, U{x = 0) =0. Here the dome width in the bulk 
is 6 = OAw while the dome in the groove has b = 0.5w 
resulting in whiskers that extend out of the bulk dome by 
25%. Experimentally, in contrast, both the FIB etched 
groove and the effective grooves along the JV stacks have 
a much stronger effect: the FIB groove in Fig. [9] shows 
similar extent of the whiskers for groove depth of just 

6 = 0.02(i, whereas the JV stacks result in whiskers that 
commonly extend all the way to the sample edges, with 
effective groove depth of the order of just S = 0.0 Id. 

The above calculation, however, does not take into ac- 
count the current Ix that flows along the groove-bulk 
interface, see Fig. [TTJ which may significantly enhance 
the described effect. This current plays a similar role to 
that of the edge current in the GB model and may sig- 
nificantly enhance the magnetic induction in the groove. 
Vortices can penetrate from the groove into the bulk only 
if the difference in their line energy cqS is canceled by the 
Lorentz potential due to that part of Ix which flows in a 
region of width ~ S near the groove, 

^ Jx{y)dy = eoS, (2) 
c Ji 

where Jx{y) is the appropriate sheet current. Jx{y) al- 
ways keeps the induction in the groove higher than in 
the bulk, thus enhancing the effect of the groove and 
increasing the extent of the vortex whiskers. In the 
Meissner state Jx{y) = and no current flows along 
the groove edges. At low fields, when the vortices are 
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Figure 10. (color online) The equilibrium potential curves 
(blue), lJ(x — 0) = 0, and induction profiles (green) for two 
regions of a superconductor having different thicknesses, d 
(continuous line) and d (dashed line), assuming both regions 
can be described by the GB mechanism. The difference in the 
edge barrier height in the two regions which is proportional 
to ^/d = {d — d)/d^ is of the same order of magnitude as the 
normalized dome width difference {h — b)/b. 
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Figure 11. A numerical solution for the sheet current stream 
lines of a two-dimensional rectangular film of half-length L 
and half- width W {W/L = 0.4) with two grooves (shaded 
regions), each 21 = 0.06L wide, at H^/Hci = 0.44. The 
arrows indicate schematically the current flow directions. It 
is assumed that the penetration of the vortices into the groove 
begins when the sheet Meissner current at a distance d/2 from 
the platelet edge reaches cHci/2ti. Here we took d = 0.05L. 



present only in the groove, Jx{y) extends over the entire 
bulk. With increasing field, dome is formed in the bulk, 
decreasing Jxiv)- When the field is high enough so Jx 
flows only in the region of width ~ (5 at the interface, 
the total interface current Ix reaches its minimum value 
= ceoS/^o = cHciS/47T given by Eq. ([2]). The differ- 
ence in the inductions at the groove center and bulk ABgr 
which is proportional to Ix , depends on the geometry of 
the groove. For example, if the whole current Ix flows 
near the groove, one has ABgr = {^T^/cjlx/d for the nar- 
row groove, / <C and ABgr = (4/c)/a^/nn the opposite 
limiting case, / ^ d. Note that the enhancement of the 
field in the groove is essentially a demagnetization effect 
due to partial expulsion of the flux from the bulk. It is 
clear that for very wide grooves the effect is negligible 
and the result of Fig. [TOl holds. 

The sheet current J{x^y) and magnetic induction 
Bz{x^y) can be calculated numerically for a finite size, 
two-dimensional film using the method described in 
Ref. 0. This method can be applied to the problem 
at hand in the limit of d <^ I <^ w. Here we demonstrate 
by this calculation that the dome in the groove can be 
noticeably wider than the dome in the bulk. To carry out 
the calculation, we rewrite condition (|2]) using the known 
solution for the sheet current and magnetic induction in 
the case when there is a magnetic flux in a narrow slit 
of the width 21 between two wide infinitely long strips 
in the Meissner state^^ This solution can be used only 
if Jx flows in the region essentially wider than i.e., if 
{4:71 /c)Ix ^ Hci^- Then, Eq. (|2j) can be rewritten as 

B,{x = 0,y = 0) = AH,i, (3) 

where A is a factor which depends on the distribution of 
the current density in the layer of the thickness S near 
the groove. Strictly speaking, this distribution is essen- 
tially three-dimensional and cannot be found from the 



solution of Ref. 'so'. We estimate that the factor A lies in 
the interval from ^/SJ^ to ^/d/L Since the calculation is 
done for a strictly 2D film, S enters the calculation only 
through A which we take as A = 0.5 in the figures shown 
next. Figures [11] and [121 show the current stream lines 
and magnetic induction of a 2D superconductor having 
two grooves of width 21 = 0.06L at H^/Hd = 0.44. The 
interface current Ix as well as the flux focusing into the 
slits is clearly seen. Figure [12] shows that the dome in 
the grooves is wider than the bulk dome forming whiskers 
protruding into the flux free region. The induction in the 
grooves is much higher than in the bulk and interestingly, 
it is even higher than the applied field. The paramagnetic 
response of the grooves, Bz > H^/is similar to the situa- 
tion described in Ref. 48, where such a phenomenon was 
described in terms of demagnetization effect. 

Our numerical calculations also provide a qualitative 
understanding of the behavior on the descending field as 
follows. When the applied field is decreased, the bulk 
dome and the whiskers along the grooves expand un- 
til the whiskers reach the sample edges. Since the bulk 
and groove vortices remain in equilibrium, on further de- 
crease of Hz the vortices are drained out from the bulk 
through the groove. This process results in reduction 
of the width of the hysteresis loop and in magnetiza- 
tion value that is closer to the equilibrium magnetization 
on the descending branch. The inductions at both the 
groove and bulk centers decrease with decreasing Hz, 
but their difference, ABgr oc I^ = cHciS/Air remains 
constant. Thus as Hz is reduced to zero the bulk dome 
contracts until all its vortices are drained out at a field 
Ho = Hij=o such that the induction at the groove center 
equals ~ ABgr. In order to estimate i^o, we note that 
at this field Bz = in the bulk and the sample is in 
the Meissner state everywhere except in the grooves. In 
the Meissner state the total screening current circulat- 
ing in a strip is /tot = cHqw /2t:^^ Near the grooves 




Figure 12. (color online) A numerical solution for magnetic 
induction Bz corresponding to the conditions of Fig. [TTl The 
field in the two grooves is higher than the applied field similar 
to the situation described in Ref i48i. 
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this current splits into two paths: flows parahel to the 
groove and /edge flows across the groove edge (see Fig. 
[TT]) . Thus, we obtain Hq ^ HciS/2w + 27r ledge /cw. We 
can thus evaluate Hq by noting that the edge current lies 
in the interval < /edge ^ cHcid/ Ai:.^ This estimate for 
the range of Hq is in reasonable agreement with the ex- 
perimental values of Hq seen on the descending branches 
ofb{H^) in Fig. [3b. 

The above description provides a qualitative explana- 
tion of the creation of vortex whiskers, the suppression 
of hysteresis, and the shrinkage of the bulk vortex dome 
on descending fleld. Quantitatively, however, the groove 
depths that are required in order to describe the experi- 
mental data are still about an order of magnitude larger 
than the ones obtained from Eq. ([T]). Our numerical 
calculations are applicable only for grooves that are not 
too narrow d <^ I <^ while the JV stacks create effec- 
tive grooves with I <^ d. In this case we expect further 
enhancement of the demagnetization effects that could 
perhaps resolve the discrepancy. Interestingly, the JV- 
PV crossing energy ex was evaluated experimentally to 
be about seven times larger than the theoretical value 
by comparing ex with the vortex line energy in wide 
grooves. The origin of this discrepancy could be re- 
lated to the same observation that the effective grooves 
due to JV's are much narrower, thus having a large fleld 
enhancement due to demagnetization effects in the thin 
samples in perpendicular fleld. 

The presented model, however, does not explain why 
the penetration fleld decreases signiflcantly with increas- 
ing Hx. In the full Meissner state Ix = and therefore 
Hp cannot be affected by it. Since in the GB model 
Hp {2Hci/7r)^/d/w, the bare reduction of the thick- 
ness along the JV effective grooves should lower Hp in 
our case by less than 1%, in sharp contrast to more than 
50% reduction in Hp in Fig. [31(b). This essential reduc- 
tion of Hp may be due to a more intriguing mechanism in 
the vicinity of the edges in which an inhomogeneous pen- 
etration of the vortices into the edge region leads to a con- 
centration of the magnetic fleld in the groove edges. As a 
result, an enhanced penetration of the vortices along the 
groove develops that noticeably decreases the observed 
penetration fleld. 

VI. SUMMARY 

We used the DM0 technique to image the fleld distri- 
bution and the local magnetization of thin BSCCO crys- 
tals at elevated temperatures at which the magnetic hys- 
teresis is governed by the geometrical barrier mechanism. 
In this regime a vortex dome is present in the central part 



of the sample surrounded by a vortex-free potential bar- 
rier region. We found that stacks of Josephson vortices 
caused by an in-plane magnetic fleld signiflcantly reduce 
the GB and allow formation of pancake vortex chains 
or whiskers in regions that should otherwise be vortex- 
free. These whiskers extend from the dome up to the 
sample edges, thus forming easy channels for vortex flow 
through the potential barrier. As a result the magnetic 
hysteresis is reduced and the magnetization loops become 
reversible at elevated in-plane flelds. We also found that 
a FIB etched narrow and shallow groove has a similar 
effect on the GB. 

Our analysis shows that the bare reduction of the vor- 
tex energy due to the JV-PV crossing energy along the 
JV stacks or due to reduced thickness along the etched 
grooves is one to two orders of magnitude lower than the 
height of the GB, and therefore cannot account for the 
observed phenomena. We present a model in which the 
effect of the grooves is signiflcantly enhanced due to the 
demagnetization factor in the platelet geometry in per- 
pendicular fleld. The flux that is partially expelled from 
the pristine regions is focused into the grooves, thus en- 
hancing the local induction in the grooves and suppress- 
ing the barrier. The presented numerical calculations 
demonstrate qualitatively this demagnetization enhance- 
ment, and provide an important insight into the long- 
standing puzzle of suppression of the macroscopic GB 
by in-plane fleld. Our numerical method can only treat 
grooves that are wide on the scale of sample thickness, 
which are much wider than in the experimental situation. 
A proper calculation of narrower grooves is therefore re- 
quired, and is expected to provide a better quantitative 
agreement with the data. One important experimental 
observation, the suppression of the penetration fleld by 
the JV stacks, however, remains unexplained within the 
model. This feature is apparently a result of microscopic 
inhomogeneities in the sample edge regions introduced 
by the JV's that facilitate vortex penetration into the 
sample. Such inhomogeneities cannot be treated within 
our current numerical studies and will be the subject of 
future investigations. 
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